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Abstract 

We report the status of a high-statistics Monte Carlo simulation of 
non-self-avoiding crystalline surfaces with extrinsic curvature on lat- 
tices of size up to 128^ nodes. We impose free boundary conditions. 
The free energy is a gaussian spring tethering potential together with 
a normal-normal bending energy. Particular emphasis is given to the 
behavior of the model in the cold phase where we measure the decay 
of the normal-normal correlation function. 
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1 Introduction 



In recent years there has been a lot of interest in the statistical mechanics of 
crystalline and fluid surfaces [Q]. The former is believed to describe physical 
polymerized membranes [|2| and the latter may be a regularization of string 
theory. 

We focus our study on crystalline surfaces with bending rigidity embed- 
ded in R^. It is conjectured that this model has a second order phase transi- 
tion driven by the competition between entropy and the bending energy [^. 
The high temperature phase is characterized by crumpled configurations. In 
the low temperature phase the system is no longer isotropic and the surfaces 
are roughly flat. 

One may wonder what stabilizes the flat phase. The theoretical argument 
is that the in-plane elastic constants prevent the surface from fluctuating 
arbitrarily in the embedding space. This leads to an effective long wave 
stiffening of the surface [Q, ^ . 

The crumpling transition has been studied numerically with simulations 
explicitly incorporating 2-d elastic constants [Q, ^, There have also been 
simulations with a simple gaussian spring potential playing the role of the 
tethering potential §, |, ||, 0, [ll|, In both cases evidence has been 
presented for a continuous phase transition. 

In the latter class of models the equilibrium spring length is taken to be 
zero, and a simple calculation indicates that the microscopic elastic constants 
vanish |jl3[. It is tempting to argue, therefore, that the flat phase of these 
models is not truly stable, even in the limit of large bending rigidity. 

Our ultimate aim is to carefully compare the behavior of the appropriate 
observables in the cold phase as a function of the equilibrium spring length. 



2 The Model 

Consider a system of particles connected to form a triangular 2-d mesh 
embedded in 3 dimensions. Let each particle be labeled by an internal 
discrete coordinate system x = (xi,X2) denoting its position on the mesh. 
Its actual position in the embedding space is given by the 3 dimensional 
vector r(a;i,X2). The action has a tethering potential and a bending energy 
term. Our choice is to use simple gaussian springs between the vertices as a 
tethering potential and a normal-normal interaction as the bending energy 
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term. Therefore the action is 

5= E (^xx')' + A^ (l-n,-n^). (1) 

(xx'> (a/3> 

Here the subscripts x, x' label the vertices and /xx' is the distance between 
the vertices x and x' in the embedding space. The subscripts a, /3 label the 
faces (triangles) of the surface, n is the unit normal to the face and A is the 
bending rigidity. The sums extend to nearest neighbours. Eq. |^ describes 
phantom surfaces since it does not include self-avoidance. The surface is a 
rhombus with free boundaries cut out of a triangular lattice. In the case of 
a spring of length a, (/xx') has to be replaced by (/xx' — a). In this case our 
model would closely resemble the one of 0] discussed above. 

We focused our analysis on the following observables: the specific heat, 

C. = y {{S!) - (Se^^^ 

Here Se is the bending energy term considered above and V is the total 
number of vertices. 

The radius of gyration, 

^^ = ^(Er^r^)• (3) 

i 

Here is the position of the node i in the embedding space referred to the 
center of mass. This observable measures the physical extent of the surface 
and its scaling behavior with system size defines the size (Flory) exponent 
u, via the relation Rg oc L'^ . The exponent u is related to the Hausdorff 
dimension (Ih via the relation z/ = 2/(1^. 

The eigenvalues of the inertia tensor; these eigenvalues give information 
on the shape of the surface and how it scales with system size. They are 
obtained by diagonalizing the anisotropic part of the inertia tensor 

Aat = Y.^'a{iK{i) (4) 

i 

where a, b refer to the components of the vector r'. 
The normal-normal correlation function, 

GiR) = {l E i^a-no)). (5) 

Here the sum extends to all triangles of the surface which have a geodesic 
distance R from the center of the surface O. The angle brackets represent 
the Monte-Carlo average. 
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Figure 1: The specific heat as a function of the bending rigidity. 

3 Theoretical Predictions 

A self-consistent perturbation theory analysis of the continuum model [^] 
yields predictions for the critical exponents. The exponents of interest are 
the size (Flory) exponent u and the roughness exponent ^. The roughness 
exponent is defined by the scaling of the minimum eigenvalue of the inertia 
tensor (^), 

A oc L^^. (6) 

At the critical point the theory predicts = ^ = 0.732 while in the 
cold phase 1^ = 1 and ^ = 0.59. 

As far as the normal-normal correlation is concerned the only analytical 
result is for A = (gaussian model). In this case the correlation function 
follows a decay law G{R) oc —1/R'^. 
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Figure 2: The scaling of the radius of gyration squared with system size in 
the cold phase(A = 1.1). The fit gives a value of = 1.05(4). 



4 Numerics and Results 

We performed Monte-Carlo simulations of systems of sizes 16^ to 128^ ver- 
tices. We used the single hit Metropolis algorithm. The largest lattice was 
simulated on a MASPAR MPl massively parallel processor, while all other 
sizes were simulated on workstations. We gathered statistics of the order of 
30-50x10^ sweeps per data point for the largest lattices (64 and 128). Our 
statistics are comparable for the smaller lattices. 

As can be seen in Fig. |l] the specific heat shows a growing peak with 
system size. Presently our statistics are not yet sufficient to allow a reliable 
estimate of the exponent a which characterizes the growth. Preliminary fits 
indicate a value of a ~ 0.5 consistent with the value obtained in |jl^, ^] 
using the same method. Our best estimate for the critical value of the 
coupling is around A ~ 0.79. Work is currently under way to gather better 
statistics and perform a Ferrenberg-Swendsen type analysis. 
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Figure 3: The normal-normal correlation function at A = 0.8 (around the 
phase transition). 

Fig. ^ shows the radius of gyration versus system size at a fixed value of 
A (1.1). The data fits well to a scaling ansatz with v = 1.05(4), as expected 
in the flat phase. In the crumpled phase the data does not fit a power law 
behavior, indicating, as expected |5|, a log-like scaling [du = oo). Our 
estimate of the critical coupling is not precise enough to allow for a fit to 
Rg at the transition. 

Figures ^ and ^ show the normal-normal correlation around and above 
the phase transition respectively. Fig. ^ demonstrates the effect of finite- 
size corrections to the value of the critical coupling. For the smaller volume 
L = 32 the correlation decays to zero with r, but fits indicate a non-zero 
asymptotic value for L = 64. One possible reason is that, due to the volume 
dependence of the pseudo-critical coupling, the smaller volume is in the 
crumpled phase while the larger is in the cold (flat) phase. Note that for 
large r the correlations will always decay to zero because of our choice of 
boundary conditions. This is a flnite size effect and the data close to the 
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Figure 4: The normal-normal correlation function at A 
phase). 



1.1 (in the cold 



boundary has to be excluded from the fits. 

Figure ^ shows the correlation in the cold phase. The data for L = 64 
fits well to a behavior 



(7) 



with a = 0.51(1) and h = 0.3(1). Data for higher values of A and of the 
system size L show a consistent behavior. 

This is very important — the presence of a non-zero asymptote for the 
normal-normal correlation function indicates that the normals remain or- 
dered on a macroscopic scale. 

This result supports the existence of a stable flat phase. Our present 
focus is on the precise nature of this phase; in particular we would like to 
know if there is a well defined roughness exponent ^ and if so how it depends 
on the bending rigidity. The existence of a flat phase with a roughness 
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exponent independent of A would indicate that the entire flat phase A > Ac 
is critical and that this model is in the same universality class as ^ . 

Another possibility is that the roughness exponent will depend on A. In 
this case the flat phase would still be stable but not critical. 

A possible interpretation of the stability of the flat phase is that the 
model discussed here dynamically generates a non-zero length scale which 
can serve to define non-zero renormalized elastic constants. This simple 
model could then be used to study the properties of physical (polymerized) 
membranes in regimes in which self-avoidance is irrelevant. 
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